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Abstract 

We  present  wellposedness  results  for  a  multi-dimensional  electromagnetic  interro¬ 
gation  problem  involving  Maxwell’s  equations  in  variational  form  coupled  with  general 
polarization  equations  and  perfectly  matched  layers  (PMLs)  at  finite  boundaries. 


1  Introduction 

In  this  paper  we  present  theoretical  foundations  (existence,  uniqueness,  continuous  depen¬ 
dence)  for  a  particular  two  dimensional  Maxwell  system.  This  system  is  a  multi-dimensional 
version  of  the  one  dimensional  systems  treated  in  [2,  Chapter  3].  The  2-D  problem,  carefully 
formulated  and  treated  numerically  in  [1],  arises  from  pulsed  microwaves  emanating  from 
a  semi-infinite  strip  of  a  finite  antenna,  producing  obliquely  incident  waves  on  a  dielectric 
material  target  as  depicted  in  Figure  1.  The  problem  is  formulated  with  distributional  (in 
space)  input  currents,  a  general  convolution  polarization  relationship  that  includes  most 
models  for  dielectrics  currently  found  in  the  literature,  and  a  finite  computational  domain 
with  partially  absorbing  boundary  conditions  as  well  as  perfectly  matched  layers  (PMLs)  to 
attenuate  incident  energy  at  the  boundaries.  The  resulting  system  involves  a  second  order 
(in  time)  Maxwell’s  equation  in  variational  form,  coupled  with  the  polarization  convolution 
and  PML  constraint  equations  most  conveniently  formulated  in  terms  of  distributions  and 
convolutions.  While  we  follow  the  general  approach  to  a  theoretical  framework  developed  in 
[2],  the  multi-dimensional  aspects  of  the  electric  field  as  well  as  that  of  the  computational 
domain  with  PML  boundaries  present  new  theoretical  challenges.  We  provide  arguments 
here  to  successfully  overcome  these  difficulties  and  establish  new  wellposedness  results. 
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2  Formulation  of  Problem 

In  this  section  we  provide  an  overview  of  the  derivation  of  the  wave  equation  of  interest  given 
in  [1],  Maxwell’s  equations  for  a  region  with  free  charge  density  p  =  0  are  given  by 

V  •  .D  =  0 
V-B  =  0 
Vx£  =  —dtB 
V  x  H  =  dtD  +  J, 

— #  — t  — t 

where  the  vectors  in  (1)  are  functions  of  position  r  —  (x.  y.  z )  and  time  t.  and  J  —  Jc  +  Js, 

— t  t 

where  Jc  is  the  conduction  current  density  and  Js  is  the  source  current  density.  We  assume 
only  free  space  can  have  a  source  current,  and  thus  either  Jc  =  0  or  Js  =  0,  depending  on 
whether  or  not  the  region  is  free  space. 

In  order  to  allow  for  a  diagonally  anisotropic  medium,  we  use  the  diagonal  operator 

T  =  diag  (Tx,Ty,Tz). 


(1) 


Here,  given  /  e  C[ 0,  oo),  T*/  =  gi  *  /  for  i  =  x,  y,  z.  where  gi  is  a  tempered  distribution  as 
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described  in  [1],  The  diagonally  anisotropic  material  constitutive  relations  are  then  given  by 


D  =  g*  (TE) 

B  =  g0  TH  (2) 

Jc  =  aTE1 


where  g  is  a  tempered  distribution,  convolution  is  in  the  time  variable,  g0  is  the  magnetic 

permeability  of  free  space,  and  the  constant  a  >  0  is  the  conductivity  arising  in  Ohm’s  law. 

By  restricting  our  attention  to  a  region  on  which  T  and  g  are  constant,  using  (2) 

— # 

and  (1),  and  defining  the  electric  polarization  vector  P  by 

g*E  =  eo  erB  +  P-  (3) 

where  er  is  the  relative  permittivity,  we  obtain  the  wave  equation 

goto +  go P  +  F o&E  —  (V  •  ( T\7))E  =  —goT  1JS,  (4) 

where  '  =  dt.  "=  df,  and 

f  =  {TxTyTzylT.  (5) 

We  now  express  equation  (4)  in  variational  form.  Let  f l  have  a  piecewise  smooth 
boundary,  and  let  (f>  be  a  test  function  on  Q.  For  a  function  /  on  O  define 

(fA)=  f  f  dV. 

Jn 

If  /  and  (f>  are  also  defined  on  let 

(fA)dn=  f  f{r)<p(r)ds. 

Jon 

Let  u  equal  one  of  the  standard  basis  vectors  x ,  ig  or  z,  and  define 

E  =  E  ■  u 

P  =  P-u  (6) 

Js  =  (T-1  J8)  •  u. 

The  variational,  or  weak,  formulation  of  (4)  is  then  given  by 

goe0(erE1  </>)  +  go(P,</>)  +  g0(aE1  </>)  +  {d^ E),  drf)  -  ^ 

i=x,y,z  i=x,y,z 

=  -go{jsA)-  (7) 
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Note  that  T  1  has  disappeared  on  the  right  side  of  (7).  This  is  because  in  our  problem  a 

current  source  will  only  exist  in  free  space  where  T  =  I.  Hence  the  right  side  of  (4)  reduces 
— # 

to  either  0  or  —/j,0Js. 

Our  computational  domain  V  will  be  the  union  of  a  finite  number  of  regions  flk  on 
which  T  and  g  are  constant  in  the  spatial  variables.  Hence,  the  components  of  a  solution 
of  Maxwell’s  equations  with  constitutive  relations  (2)  will  satisfy  (7)  on  each  £lk.  To  extend 
to  the  entire  domain  we  must  first  determine  V.  As  in  [1],  we  assume  that  E,  P,  and  Js 
are  independent  of  y.  so  our  computational  domain  is  contained  in  the  xz-plane.  Let  X 
and  Z  be  closed  intervals  in  the  respective  x-  and  2-axes,  and  let  V  —  X  x  Z  denote  the 
computational  domain. 

We  partition  the  interval  X  into  disjoint  closed  intervals  X  .  X0,  and  X+  such  that 

max  X  =  min  A'0 
max  W  =  min  A+, 

and  partition  the  interval  Z  into  disjoint  closed  intervals  Z _,  Zv.  and  Zd  such  that 

max  Z  =  min  Zv 
max  Zv  =  min  Zd . 

The  regions  Zv  and  Zd  are  the  vacuum  and  dielectric  regions  respectively.  Let 


Zq  —  Zv  U 


The  region  X0  x  Zq  will  be  our  domain  of  interest.  The  buffer  region  V  \  (A0  x  Zq)  outside 
our  domain  of  interest  will  contain  the  PMLs.  This  problem  geometry  is  illustrated  in  Figure 
2. 

We  also  need  to  choose  an  operator  T.  In  order  to  guarantee  zero  reflection  at  the 
PML  interfaces,  T  must  be  chosen  carefully.  In  [1],  T  is  taken  to  be 

(S^s,  0  0  \ 

T  =  0  SXSZ  0  ,  (8) 

v  0  0  sxsz  7 

where  Sx  and  Sz  are  convolution  operators  on  C[ 0,  oo)  with  kernels 

ux(t,  x )  =  6(t)  +  Pxixfflty), 

Uz(t,  z)  =  5(t )  +  pz(z)n{t). 

Here  5  is  the  Dirac  delta  function,  H.  is  the  Heaviside  function,  and  both  /3X  and  /3Z  are 
piecewise  constant,  convex  functions  such  that  /3X  =  0  in  X0  and  f3z  =  0  in  Zq.  It  is  a 
straight  forward  exercise  to  show  that 

f  =  diag(Sx2,/,S;2),  (10) 
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Figure  2:  Computational  domain. 


where  /  is  the  identity  operator,  and  Sx  1  and  Sz  1  are  convolution  operators  with  kernels 

ux\t,  X )  =  S(t)  -  A c(x)/H{t)e^x)t,: 

u^it.z)  =  5(t )  -  /4(z)^(*)e~^(z)f. 

For  computational  reasons,  14  and  Vy  are  introduced  and  are  given  by 

Vx  +  E  =  TXE 
Vz  +  E  =  TZE. 

Then,  given  the  above  comments,  in  any  region  Q  in  which  j3x  and  /3Z  are  constant,  the 
variational  equation  (7)  can  be  written 

Hoe0{erE,  4>)  +  /4>(A  4>)  +  l^o(o'E1  cf>)  +  {dx(E  +  14),  dx<j) )  +  {dz(E  +  14),  dz(f> ) 

~  {dx(E  +  14 ),nx<p)dn  -  (dz(E  +  14),  nz(p)dn  =  -/j0(  0).  (13) 

Provided  that  the  PML  parameters  /3X  and  (3Z  are  chosen  carefully,  it  is  argued  in  [1] 
that  (13)  will  hold,  to  a  sufficient  degree  of  accuracy  for  computational  purposes,  on  the 
computational  domain  V.  Then,  the  variational  equation  that  is  used  for  computations  is 
given  by 


Hoe0(erE,  <p)  +  <t>)  +  Ato(cr-E',  <t>)  +  {dx(E  +  14),  dx<p)  +  ( dz(E  +  14),  dz<fi) 

~  (dx{E  +  14),  nx(j))dv  -  (dz(E  +  14),  nz<\>)dv  =  -//0{  Js,  </>).  (14) 
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Finally,  we  add  absorbing  boundary  conditions.  Let  C  Xl  C+x.  C-Zi  and  C+z  de¬ 
note  the  four  boundaries  of  V  with  outward  normals  in  the  —x,  +x,  —z  and  +z  directions 
respectively.  We  use  the  boundary  conditions 

(/-4  +  dt)E  =  cdxE  on  C  x 
(j%  +  dt)E  =  —cdxE  on  C+x 
{fiz  +  dt)E  =  cdzE  on  C  z 
E  =  0  on  C+z 

where  c  =  l/yAk^o-  The  boundary  condition  on  C+z  models  a  dielectric  with  a  supracon- 
ductive  backing.  It  is  shown  in  [1]  that  it  is  reasonable  to  expect  that  each  of  the  remaining 
three  boundary  conditions  will  reflect  only  approximately  3%  of  incoming  energy.  Equation 
(14),  with  a  rescaling  of  the  time  variable  by  a  factor  of  c  =  l/y'eo/io  and  a  rescaling  of  P 
by  a  factor  of  1/eo,  then  becomes 


{erE,  <j>)  +  (P ,  <J>)  +  rjo(aE,  (f>)  +  (dx(E  +  Vx),  dx(f>)  +  (dz(E  +  Vz).  dz4> ) 

+  ((“Ax  +  dt)(E  +  14),  4>)c-xuc+x  +  ((-As  +  dt)(E  +  Vz),  <p)c-z  —  ~Vo (Ts,  <f>),  (16) 
where  7]o  =  \/Ao  Ao  is  the  impedance  of  free  space. 


3  Wellposedness 

In  this  section,  we  prove  that  the  initial  value  problem  for  the  variational  equation  (16)  is 
wellposed.  That  is,  given  initial  conditions 

E(Q,x,z)  =  3>(x,z),  E(0,  x,  z)  =  \k(x,  z)  (17) 

we  will  show  that  there  exists  a  unique  solution  of  (16),  (17)  that  depends  continuously  on 
$,  4/  and  Js.  Without  loss  of  generality  we  take  er  =  1  in  (16). 

We  now  explicitly  define  the  space  V  of  test  functions  <f>.  This  space  is  given  by 
V  =f  HlR(V )  =  {(f)  ^  Hl(V)  |  (f>(x,z )  =  0  whenever  ( x,z )  G  C+z},  where  if1(X>)  =  {(f)  G 
L2(V )  |  V  (f)  G  L2(V)}.  Following  the  approach  in  [2],  we  seek  solutions  t  — »  E(t)  with 
E(t)  G  V ,  where  t  G  [0,  T],  We  do  this  in  the  context  of  a  Gelfand  triple  V  •—>  H  V*  (see 
[7]).  Such  an  E  is  known  as  a  weak  solution  of  (16)  and  is  required  to  satisfy  E  G  L2(0,  T;  H ), 
E  G  L2(0,  T ;  H)1  and  E  G  L2(0,T;  V*)  =  L2(0,  T;  G)*.  We  therefore  interpret  _(•,•)  to  be 
(•,  -)h  in  all  of  the  terms  of  (16)  except  for  (E(t),  <f>(t))  and  (J{t),(f){t))  =  (— r7o«/s(f),  </>(i)}, 
in  which  case  (•,  •)  is  given  by  the  duality  product  (•,  • )v*,v ■  We  allow  J{t)  values  in  V* 
since  in  our  application  (see  [1])  it  is  of  the  form  g(x,t)S(z). 

To  rewrite  (16)  in  the  form  with  which  we  will  work,  we  first  express  Vx  and  Vz  in 
convolution  form,  namely 

v,(t)  =  H„E(t)  =  g,(t)*E(t),  (lg, 

V  (t)  =  H„E(t)  =  gz(t)*E(t), 
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where 


(19) 


9x(t)  = 

9z(t)  = 

The  convolution  kernels  uxl  and  u^1  are  defined  in  equation  (11). 

For  notational  convenience,  we  define 

dlcig(HgxJ  Hgz  ). 

Then,  given  that  Hgx  and  Hgz  are  constant  in  x  and  z  respectively  on  the  subregions  of 
V  mentioned  above,  from  (18)  we  have 

(dxVx,  dx<j>)  +  {dzVz,  ds<t>)  =  (HXjZV  E,  V  </>).  (20) 

By  definition,  /3X  is  constant  on  C  x  and  C+x ,  and  j3z  is  constant  on  C_z.  In  our 

computations,  the  PML  region  is  constructed  so  that  the  values  of  /3X  on  C  x  and  C+x  and 

•  •  def 

j3z  on  C  z  are  the  same.  We  will  refer  to  this  value  on  the  boundary  C  =  C  x  U  C+x  U  C  z 
by  j3c.  Our  assumption  on  /3X  and  /3Z  also  has  the  result  that  the  operator  Hgx  on  C  x  and 
C+x  is  identical  to  Hgz  on  C  In  the  sequel,  we  will  denote  this  operator  on  C  with  kernel 
9c  by  Hgc. 

Using  Liebniz’s  rule  and  the  fact  that 

P(t,x,z)=  /  g(t  —  s,  x,  z)E(s,  x,  z)ds  (21) 

Jo 

we  obtain 

P(t,x,z)  =  /  g(t  —  s,  x,  z)E(s,  x,  z)ds  +  g(0,  x,  z)E(t1  x,  z)  +  <?(0,  x,  z)E(t1  x,  z).  (22) 

Jo 

We  assume  that  g,  g ,  and  g  are  bounded  on  [0,T]  x  V. 

Then,  given  X(x,  z)  =f  Iv(%,  z)g( 0,  x.  z),  where  Tp  is  the  indicator  function  on  V,  there 
exists  k  >  0  such  that  A  =  A  +  A;>e>0.  We  define  a  sesquilinear  form 

0i(<P,  ip)  =  (Mi  ip)  +  (V  <p,  V  ip)  for  <p,ipeV.  (23) 


It  is  readily  seen  that  is  U-continuous  and  U-elliptic.  That  is,  there  are  positive  constants 
Ci  and  c2  such  that 

01  <  Cz\<P\v\*P\v  for  all  <p,ipeV, 
oMi  <P)  >  Cl  101  v  for  (j)  e  V. 

Finally,  combining  (23),  (22),  (20),  (18),  and  (16)  we  obtain 


(E(t),  <p)  +  (7 E(t),  <p)  +  <ti (.E(t),  <p)  +  (  f  a(t  —  s)E(s)ds,  <p)  +  ( 

Jo 


H XyZV  E(t),V  <j>) 

Jo 

+  ^(H„E(t),4,)c  +  ^(E(t),  0>c  +  {H„E(t),  6)c  +  (E(t),  cf)c  =  <  J(t)A )  +  (kE{t),  6), 


(25) 
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where 


a(t,x,z )  =  Iv(x,z)g(t,x,z )  and  l{x,z)  =  Iv(x,  z)(rj0a(x,  z)  +  g(0,x,z)). 

Since  g  and  g  are  bounded  on  [0,T]  x  V  and  a  is  assumed  bounded  on  V,  so  are  a(t )  and 
7.  We  now  prove  wellposedness  for  (25),  (17),  which  implies  wellposedness  for  (16),  (17). 

Existence 

In  order  to  prove  the  existence  of  solutions  for  (25),  (17),  where  <3?  €  V  and  T  G  77,  we 
follow  closely  the  general  framework  found  in  [2,  4,  5].  We  begin  by  choosing  a  linearly 
independent  subset  {wi}^  that  spans  V,  and  we  let  Vm  =f  span{w1; . . . ,  wm}.  Define 
Galerkin  “approximates” 

m 

Em(t,x,z )  =  y ^e™(t)wi(x,z).  (26) 

i= 1 

The  coefficients  {e™}  are  determined  by  solving  the  linear  system  that  results  from  replacing 
E  by  Em  and  (j)  by  Wi  for  i  =  1, . .  .^rn  in  (25),  and  taking  Em( 0)  =  <3>m,  Em( 0)  =  \km,  where 
<3? m,  G  Vm  such  that  <3?m  — >-  <3?  in  V,  — >-  T  in  H.  Taking  <p  =  Em  in  equation  (25)  we 

then  obtain 

~\Em{t)\2H  +  \^Em(t)\2H  +  \^°i(Em(t),  Em(t))  +  (Jo  a(t  -  s)Em(s)ds,Em(t)) 

+  ^(HX)ZV  Em(t),\7  Em(t ))  —  {(^(H^V  Em(t)),  V  Em(t ))  +  —  —  (EJgcEm(t),Em(t))c 

-  ^(dt(HgEm(t)),  Em(t))c  +  jt(HgcEm(t),  Em(t))c  -  (dt(HgEm(t)),  Em(t))c 

+  2^df\Em{t)\L2(c)  +  \Em(E)  Il2(c)  =  (77(f), -Em (f))  +  (kEm{t),  Em(t)).  (27) 

The  convolution  kernel  gc  is  given  by 

9c(t)  =  -2  pcU{t)e-Pct  +  P2cm(t)e-M.  (28) 

The  kernels  gx  and  gz  are  given  by  (28)  with  /3C  replaced  by  j3x{x)  and  Pz(z)  respectively. 
Using  (28)  and  integration  by  parts,  we  obtain 

(■ HgcEm(t ),  Em(t))c  =  {~Pc9c(t)  *  Em(t),  Em{t))c  +  Pl{^H{t)e~Pct)  *  Em(t ),  Em(t))c 

-  2 pc{Em(t),  Em(t))c  -  (gc(t)Em( 0),  Em{t))c 

=  (~Pc9c(t)  *  Em(t),  Em(t))c  +  PcdW (t)e_Ajt)  *  Em(t ),  Em(t))c  ^ 

-  4 pc  f  ( Em(s ),  Em(s))cds  -  2/3c(Em(0)1  Em{ 0))c 

Jo 

~  {gc(t)Em(0),  Em(t))C. 
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Integrating  (27),  using  equations  (24)  and  (29),  and  noting  that  H^V  Em( 0)  =  HgcEm{ 0)  = 
HgcEm( 0)  =  0,  we  have 


\Em(t)\2H  +  2  [  WlEm(s)\2Hds  +  ci\Em(t)\y  +  —\Em(t)\2L2^  +  2  f  \Em{s)\2L2^ds 
Jo  c  Jo 

<  \Em{Q)\2H  +  C2|-E'm(0)|y  +  (^  +  4/3c)  | (0)  |  —  2(HXjZV  Em(t ),V  Em(t )) 

+  2ft  (5c(t)^£m(t),Em(t))c-2/5c2((?7(t)e-^)^i?m(t),Em(t))c 

+  2(</c(i).Elm(0),  Em(t))c  +  2  f  Fm(s)ds,  (30) 

Jo 

where 

^Jn(t)  =  (*7(i),  4W)  +  ( kEm(t ),  £„,(<))  +  (~  [  a(t  -  s)Em(s)ds,  Em(t )) 

Jo 

+^(dt(HgeEm(t)),  Em(t))c  +  <<9t(H^V  £m(t)),  V  Em(t)) 

T  (<9^ (Hgc Em (f ) ) ,  Em(t))c  +  8 Pc(Em(t),  Em(t))c 
=  F'rnA  (s)  +  Fm^{s)  +  Fm^(s)  +  Fm^(s)  +  Fm^(s)  +  Fmfi(s)  +  Fmj(s). 

We  want  an  upper  bound  for  the  left  side  of  (30).  To  this  end,  we  seek  upper 
bounds  for  ((H^V  £m(t),  V  £m(t))|,  |(^c(*)^m(*),  Em(t))c,  \ {(H(t)e-Pct)*Em(t),Em(t))c\, 
|  f*  Fm(s)ds\,  and  |(^c(f)Era(0),  Era(t))c|. 

To  obtain  an  upper  bound  for  |(HXjZV  Em(t),V  Em(t))\,  we  first  note  that  by  equa¬ 
tions  (18)  and  (20), 


(HX}ZV  Em(t),V  Em(t ))  =  (HgxdxEm(t),dxEm(t))  +  (HgzdzEm(t),  dzEm(t)). 
We  then  obtain  an  upper  bound  for  \(HgxdxEm(t),  dxEm(t))\. 

\(HgxdxEm(t),dxEm(t))\  <  |(  [  gx(t  -  s)dxEm(s),  dxEm(t))\ 

Jo 

—  I  /  9x(t  s)(dxEm(s),  dxEm(t))ds\ 

Jo 

[  \{dxEm(s),dxEm(t))\ds 
Jo 

[  \dxEm(s)\Hds\dxEm(t)\H 
Jo 


<  Cg 


<  Cg 


<^(  /  \dxEm(s)\Hds)  +  e\dxEm(t)\2H 
Tc 


< 


Tc2  F 

J  \dxEm(s)\2Hds  +e\dxEm(t)\2Hl 


(31) 
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where  cg  is  given  by 


Cg  =  2(3c  +  f32T,  (32) 

which  is  an  upper  bound  for  gc,  gx.  and  gz.  A  computation  equivalent  to  (31)  for  the  z  term, 
combined  with  (31)  yields 


Tc2  rl 

\(UXyZV  I  |V  Em(s)\2Hds  +  e\V  Em(t)\2H 

Tc 2  rl 

l  \Em{s)\2vds  +e\Em(t)\2v. 

The  second  inequality  follows  from  the  fact  that  |V  <j)\n  <  |^|y  for  <f>  EV. 

An  argument  analogous  to  (31)  also  yields 

Tc2 

\{gc{t)*Em(t),Em(t))c  <  ^-  J  \Em(s)\2L2{c)ds  +  e\Em(t)\2L2{cy 

and 

\(('H(t)e-Pct)-kEm(t),Em(t))c\  <  \Em(s)\2L2{c)ds  +  e\Em(t)\2L2{cy 

Furthermore,  we  have 


(33) 


(34) 

(35) 


\(dc(t)Em(0),  Em(t))c\  <  Cg\Em(0)\L2(c)  •  \Em(t)  |l2(C) 

c2 

—  ^ \Em(0)  ll2(c)  +  e\Em{t)  Il2(c)- 

We  obtain  an  upper  bound  for  |  Fm(s)ds |  by  determining  bounds  for 

I  Fmj(s)ds\,  i  —  1, . . . ,  7. 

Jo 

By  arguments  found  in  [2,  p.  30], 

\  f*  Fmjl(s)ds\  <  ±-e\J(t)\2v,  +e\Em(t)\2v+1-\jm2y,  +\\Emml 

+  Jo  |  +  2^Em^S^ 


(36) 


|>  ds, 


(37) 


J  \FmM\ds  <  !^k2\Em(s)\2H  +  ^\Em(s)\2H^ds, 

f*\FmM\ ds  <  K,  f  \Em{s)\2vds  +  K2  f  \Em{s)\ 

Jo  Jo  Jo 


2Hds. 


(38) 

(39) 
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Obtaining  bounds  for  the  additional  terms  is  tedious,  but  straightforward.  For 
Jo  |Fjti,4  (s)  \  ds  we  first  note  that 

( dt(HgcEm(t, )),  Em(t))c  =  (— 2/3c£'m(t),  Em(t))c  +  (hc(t)  *  Em{t ),  Em(t,))c,  (40) 

where 

fic(f)  =  3^2cn{t)e-^  -  f33cm{t)e-^.  (41) 

Hence, 


|(at(if9c^m(f)),^m(f))c|  <  2/3c|£'ra(t)|^2(c) 

+  |<(3 *  Em(t),  Em(t))c\, 

<  2(5c\Em(t)\L2^ 

T(3f32  +  /33T)2  rt 


(42) 


+ 


J  \Em{s)\L2(C)dS  +  2 \Em(t)  Il2(C)- 


The  second  inequality  in  (42)  follows  from  an  argument  analogous  to  (31),  with  e  =  |.  Then 
f ;  |(as(if9cifm(S)),  Era(s))c|ds  <  (4&  +  1  +  T2(3f32  +  ftT)2)/ 2  jT  |^ra(s)|i2(c)ds.  (43) 

For  ft  \Frriy5(s)\ds1  we  note  that,  analogous  to  (40), 


(dt(HXjZV  Em(t)),  V  £m(t))  =  ~(2/3xEm(t),Em(t))  +  (hx{t)  *  dxEm(t),  dxEm(t))  (44) 

-<2  pzEm(t),Em(t))  +  ( hz(t)*dzEm(t),dzEm(t )), 

where  and  are  defined  by  (41)  with  the  c  replace  by  x  and  z  respectively.  Repeating 
the  arguments  found  in  (42)  for  the  absolute  value  of  both  the  x  and  the  z  terms  in  the  right 
side  of  equation  (45),  and  noting  that  /3X,  (3Z  <  /3C.  we  obtain 

t  |(3s(HXjZV  Em(s)),  V  Em(s))\ds  <  (4/3c  +  1  +  T2(3f32  +  &3T)2)/ 2  f  |V  Em(s)\2Hds 
Jo  Jo 

<  (4/3c  +  1  +  T2(3j32  +  f33T)2)/2  f  \Em(s)\2vds. 

Jo 

t  (45) 

For  fQ  \Fmfi(s)\ds,  we  use  a  modification  of  the  arguments  found  in  (42)  to  obtain 

\{dt(HgcEm(t)),  Em(t))c\  <  2 /3c(Em(t),Em(t))c 

+  eT(3f32  +  f33T)2  J  \Em(s)\2L2^ds  +  -^\Em(t)\2L2^cy 

(46) 

Then 

j\(dt{HgcEm{s)),Em{s))c\ds  <  eK3J*  \Em(s)\2L2{c)ds  +  J*  \Em(s)\2L2{c)ds,  (47) 


11 


where  K:i  =  2f3c  +  T2( 3/32  +  /33T)2. 

An  argument  analogous  to  (36)  yields 

J  \(Em(s),Em(s))c\ds<eJ  \Em(s)\\2{c)ds  +  ^  J  \Em(s)\2L2{c)ds ,  (48) 

which  gives  us  a  bound  for  |FTO;7(s)|ds. 

Finally,  combining  equation  (30)  with  the  bounds  (33),  (34),  (35),  (36),  (37),  (38), 
(39),  (43),  (45),  (47),  and  (48),  we  obtain 


\Em{t)\ii  +  (ci  -  4e)|£'ra(f)|y  +  I  ^  -  eCi  1 


+  2 


\y^Em(s)\2Hds  +  (2  -  eC2)  \Em(s)\2L2{c)ds 

f  {Cz\Em{s)\2H  +  C^\Em(s)\y  +  C5\Em(s)\2L2(CAds,  (49) 
Jo 


A  Cm  + 


where  Ci, . . . ,  C5  are  positive  constants,  and  Cm  is  the  constant  given  by 
Cm 


Emm2H  +  (C2  +  l)|£’m(0)|2,  +  f  2 &  — - 1-  -^7  )  l-^m(O)  Il2(C) 


+ 


JWIP  + 1 JWIP  +  J‘ \j(,)\U°- 


(50) 


We  choose  e  >  0  so  that  the  left  side  of  (49)  has  only  positive  coefficients. 

Recall  that  Em( 0)  =  d>m,  and  Em( 0)  =  \km.  Hence  the  sequences  {Em{ 0)}  and 
(£’m(0)}  are  bounded  in  V  and  H  respectively  since  they  converge  to  <f>  and  T  respec¬ 
tively.  Furthermore,  we  assume  that  J  G  771(0,T;  V*).  Then,  ignoring  the  fourth  and  fifth 
terms  on  the  left  side  of  (49)  and  using  Gronwall’s  inequality  [7],  we  conclude  that  {Em}  is 
C(0,T;H )  bounded,  {Em}  is  C(0,T;V)  bounded,  and  {Em\c}  is  <7(0,  T;  E2{C))  bounded. 
Finally,  substituting  these  bounds  into  (49),  we  find  that  {Em\c}  is  T2( 0,  T ;  E2(C))  bounded. 
Hence,  extracting  subsequences  if  necessary,  we  have  that  there  exists  E  G  L2(0,T;R), 
E  G  L2(0,  T;  if),  Ec  G  L2(0,  T;  L2(C)),  and  Ec  G  L2(0,  T;  E2(C))  such  that 


Em  — >  E  weakly  in  L2(0,  T;  V), 

Em  — >  E  weakly  in  L2(0,  T;  if), 

#m|c  ->  weakly  in  L2(0,  T;  L2(<7)), 

Em\c  ->•  weakly  in  L2(0,  T;  L2(C)), 


where  F’jc  denotes  the  function  E  restricted  to  the  boundary  C. 

We  take  E  as  a  candidate  for  a  solution  to  (25),  (17).  We  first  need  to  show  that 
E  =  E,  E\c  =  Ec  on  C.  and  F’jc  =  Ec  on  C.  To  this  end  we  note  that 

Emit)  =  Em{ 0)  +  [  Em(s)ds  (51) 

Jo 


12 


and 


(52) 


Em{t)\c  —  Em(0)\c  +  f  Em(s)\cds. 

Jo 


Passing  to  the  limit  in  the  weak  H  sense  yields 


E(t)  =  U(0)  +  [  E(s)ds, 

Jo 

Ec(t )  =  Ec( 0)  +  [  Ec{s)di 
Jo 


From  (53)  we  conclude  that  E  =  E,  while  (54)  implies  that  Ec  =  Ec  on  C .  Thus,  if  we 
can  show  that  E\c  =  Ec  we  are  finished.  Noting  the  topological  equivalence  (see  [6,  p.  27]) 
between  the  standard  V  inner  product  and  the  modified  inner  product  on  V  defined  by 

(</>,  0)  d=  (01  c,  ip\c)c  +  <V  0,  V  Ip)v,  (55) 


we  see  immediately  that  weak  convergence  in  V  implies  weak  convergence  in  L2(C).  Thus 
Em  — >  E  weakly  in  L2(0,T;  V )  implies  Em\c  — >  E\c  weakly  in  L2(0,T;  L2(C)),  from  which 
E\c  =  Ec  follows  immediately. 

Using  Lemma  5.1(b)  of  [3]  we  find  that  in  fact  Em  converges  weakly  to  E  in  C(0,  T;  H), 
so  that  E  e  C(0,  T;  H)  n  L2(0,  T;  V). 

Our  candidate  E  then  satisfies 


Em  — >  E 

weakly  in  L2(0,T;  U), 

(56) 

Em  — >•  E 

weakly  in  L2(0,T]  H), 

(57) 

Em\c  — >  E\c 

weakly  in  L2(0,T;  L2(C)), 

(58) 

Em\c  — >  E\c 

weakly  in  L2(0,T;  E2[C)). 

(59) 

We  now  must  show  that  E  satisfies  (25),  (17).  We  follow  the  argument  of  [4]  and  take 
0  €  ,T]  with  0(T)  =  0,  defining  0,(t)  =f  0( t)wj1  where  {wj}^  is  defined  as  above. 

For  fixed  j  we  have  that  for  all  m  >  j,  Ern  must  satisfy 

J  { {Em (t) , ipj (t) )  +  ('yEm(t),ipj(t))  +  ai(Em(t),^j(t)) 

+  {  f  a(t  —  s)Em(s)ds,  0,(t))  +  (H^V  Em(t),  V  0,(t))  +  -(HgcEm(t),  ipj(t))c 
Jo  c 

+  ^(Em(t),'tpj(t))c  +  ^{HgcEm(t),ipj(t))c  +  ^{Em(t),ipj(t))c}dt 
=  J  +  (kEm(t),'ipj(t))}dt. 
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+ 

+ 


(60) 


Integrating  by  parts  in  the  first  term,  taking  the  limit  as  m  — >  oo,  and  the  convergences  of 
(56)-(59)  yield 

+  ('yE(t),$j(t))  +  ai(E(t),*l>j(t)) 

{  f  a(t  -  s)E(s)ds,  +  <Hx,2V  E(t),  V  iP3(t))  +  -(HgE(t),  ^(t))c 
Jo  c 

P(E(t),il’j(t))c  +  ^(HgcE(t)i‘ipj(t))c+  ^(E(t)ripj(t))c]dt 
=  +  (kE(t),ipj(t))}dt  +  0)). 

Further  restricting  ip  so  that  ip  E  C^°(0,T),  we  obtain 

f  ip(t)(-E(t),Wj)dt+ 

Jo 

J  ip(t)^E(t),Wj)  +a1(E(t),wj)  +  (  j  apt  -  s)E(s)ds,Wj)  +  E(t),W  wj 

+  (Hgc E{t),  wj ) c  +  Wj)c  +  -{HgcE(t),Wj)c  +  ~{E(t),Wj)c 


{Jiwj)  ~  {kE{t),Wj)  \dt  =  0  (61) 


for  each  Wj.  Then,  using 

jtip(t)(E(t),Wj)  =  ip(t)(E(t),  Wj)  +  ip(t)-^(E(t)1wj) 
equation  (61)  implies  that 

d  .  rl 

—  (E(t),wj)  +  (^E(t),wj)+a1{E(t),wj)  +  (  J  a(t- s)E(s)ds,Wj) 

+  (H Xi*V  E(t),  V  Wj)  +  HgcE(t ),  Wj)c  +  E(t ),  w3)c 

+  (H9fE(t).,  Wj)c  +  ( E(t),Wj)c  =  { J ,  Wj)  +  (kEm(t),Wj)  (62) 

for  each  Wj.  Since  {wj}(j)L1  spans  V  we  therefore  have  that  £7  E  L2(0,  T;  V*)  and  that  E 
satisfies  (25)  for  all  p  E  V. 

From  (51)  we  have  £7(0)  =  <f>.  To  show  that  £7(0)  =  T,  we  return  to  equation  (60). 
Integration  by  parts  in  the  first  term  together  with  equation  (62)  yields 

(E(t),Wj)ip( 0)  =  (ty,Wj)ip( 0)  for  all  j. 

It  follows  that  £7(0)  =  T,  and  hence,  £7  is  a  solution  of  (25),  (17). 
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Continuous  Dependence 

The  continuous  dependence  of  solutions  of  (25),  (17)  on  <f>,  T,  and  J  follows  from  equation 
(49).  We  set 

Km  =  \Emm2H  +  (C2  +  l)!^m(0)|y  +  |^(0)||2(C7) 

1  c1  ■ 

+  2e  I  ^l^°°(0’  TT*)  +  l^(°) \v*  +  J  \j{s)\v*ds. 

<  \Em(0)\2H  +  /Ci|Era(0)|^  +  JC2\J\2Hl{0T.v»y  (63) 

Compare  Km  with  (50).  The  inequality  (63)  follows  from  an  application  of  the  norm  rela¬ 
tionship  resulting  from  (55).  That  is,  we  know  there  exists  a  c  >  0  such  that 

\<t>\c\ l2(c)  <  cHv  (64) 

for  all  (j)  G  V.  Another  application  of  (64),  together  with  equation  (49)  yields 

\Em{t)\2H +  \Em(t)\y  <  uKm+  f  u{\Em(s)\2H  +  \Em(s)\y}ds,  (65) 

Jo 

for  some  v  >  0.  Using  Gronwall’s  inequality  again,  we  obtain 

\Em(t)\2H  +  \Em(t)\y  <  vKmevT  for  t  G  [0,T],  (66) 

We  note  that  limm_s.00  Km  <  K  where 

K  =  1*13,  +  Kill’ll,  +  K2|  (67) 

Then,  via  weak  lower  semicontinuity  of  norms,  the  convergences  (56)  and  (57),  and  (66)  we 
have 

l-®,|i2(o,T;i?)  +  I-^1l2(o,T;V)  —  vKTevT .  (68) 

Since  the  mapping  ($,  T,  J)  — »  (E,  E )  is  linear  from  V  x  H  x  ff1( 0,  T;  V*)  to  L2(0,  T;  V)  x 
L2(0,T;  H),  equations  (68)  and  (67)  give  us  the  continuous  dependence  of  solutions  (E,E) 
of  (25),  (17)  on  the  initial  data  (<f>,  T)  and  J . 

Uniqueness 

For  uniqueness  of  solutions  to  (25),  (17)  it  suffices  to  show  that  the  only  solution  of  (25) 
corresponding  to  zero  initial  data  (<f>  =  T  =  0)  and  zero  input  {J  =  0)  is  the  trivial  solution. 
Let  E  be  the  solution  corresponding  to  <f>  =  T  =  J  =  0.  Then,  passing  to  the  limit  in  (65) 
yields 

m)\l  +  I E{t)\l  <  f  H\E(s)\h  +  \E(s)\2y}ds. 

Jo 
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Gronwall’s  inequality  then  implies 


\m\2H  +  m)$  =  o. 

Hence  Eft )  e  0  on  (0 ,  T)  and  solutions  of  (25),  (17)  are  unique. 

We  summarize  the  results  of  this  section  in  the  following  theorem.  We  note  that  if  g,  g, 
and  g  are  in  L°°([0,T]  x  V),  and  a  E  L°°(T>)  as  assumed  above,  then  7  and  a  satisfy  the 
conditions  of  this  theorem. 

Theorem  3.1  Suppose  that  J  E  H1( 0,T;  V*),  7  E  L°°(V),  and  a  E  L°°(0,T;  L°°(V ))  with 
ol  and  7  vanishing  outside  ofV.  In  addition,  suppose  that  the  PML  parameters  (3X  and  j3z 
are  piecewise  constant,  convex  with  maximum  value  j3c  on  the  boundaries  C  x  U  C+x  and 
C  z  respectively.  Then  for  E  V  =  HfifD),  T  E  H  =  L2(T>),  we  have  that  solutions 
to  (25),  (17)  exist  and  are  unique.  These  solutions  satisfy  E  E  L2(0,T;H)  n  C(0,T;  H), 
E  E  L2(0, T;  H),  and  E  E  L2(0, T;  V*),  and  depend  continuously  on  (<f>,  T,  J)  as  maps  from 
V  x  77  x  if1  (0,  T;  V*)  to  L2(0,  T;  V)  x  L2(0,  T;  H). 

4  Conclusions 

The  theoretical  results  established  here  provide  a  rigorous  foundation  for  the  system  inves¬ 
tigated  computationally  in  [1],  New  arguments  to  handle  the  PMLs  in  the  context  of  a 
multidimensional  system  with  finite  antenna,  current  generated  microwave  pulses  are  given. 
Similar  formulations  should  provide  a  framework  for  a  convergence  and  error  analysis  of  the 
finite  element,  PML  computational  methods  developed  in  [1], 
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